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We have performed spin and charge transport measurements in dual gated high mobility bilayer 
graphene encapsulated in hexagonal boron nitride. Our results show spin relaxation lengths As up 
to 13 pm at room temperature with relaxation times Ts of 2.5 ns. At 4 K, the diffusion coefficient 
rises up to 0.52 m^/s, a value 5 times higher than the best achieved for graphene spin valves up to 
date. As a consequence, As rises up to 24 pm with Ts as high as 2.9 ns. We characterized 3 different 
samples and observed that the spin relaxation times increase with the device length. We explain our 
results using a model that accounts for the spin relaxation induced by the non-encapsulated outer 
regions. 
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The ability to transport spin currents over long dis¬ 
tances is a major requirement to achieve new function¬ 
alities using spin information [1]. For this purpose, 
graphene and its multilayer forms are predicted to be 
ideal materials. The low spin orbit coupling provided by 
the low atomic mass of the carbon atoms [5] together with 
their outstanding electronic properties [3] make them 
great candidates for spin transport devices. 

Even though theoretical predictions suggest spin re¬ 
laxation times (ts) up to 100 ns in single layer pris¬ 
tine graphene [4], the first experimental results for 
graphene on Si02 showed Tg = 150 ps and spin relaxation 
lengths (As) of about 1.5 pm [3, opening the debate on 
which relaxation mechanism rules the spin relaxation in 
graphene [2 IS]. Following the first experiments on sin¬ 
gle layer graphene, other experiments 01 focussed on 
the spin transport properties of bilayer graphene [9Hn]. 
These experiments reported spin relaxation times up to 
6.2 ns at low temperatures [7] and 2 ns at room temper¬ 
ature [!]■ 

A recent experiment using hexagonal boron nitride 
(BN) to encapsulate single layer graphene achieved spin 
relaxation lengths up to 12 pm at room temperature 
[12j . Another experiment, using partially suspended mul¬ 
tilayer graphene that was covered by a BN flake, achieved 
room temperature up to 3.9 ns in trilayer graphene |13j 
showing the potential of graphene/BN heterostructures 
for spin transport. 

In this rapid communication we report spin transport 
in high mobility bilayer graphene (BLG). Our samples 
consist of BLG that is partially encapsulated between 
two BN flakes. The stacks are fabricated using a dry 
transfer technique that allows us to get very clean in¬ 
terfaces without extra cleaning steps [iiiin]. After this 
process, the (BN/BLG/BN) stack remains covered by a 
PC film that is removed by keeping it in chloroform at 
50 °C for 5 hours. This step is followed by an annealing 
at 250 °C for 24 hours with an Ar/H 2 flow. The contacts 
are defined using e-beam lithography and designed with 


different widths (from 0.1 to 0.5 pm) to ensure different 
coercive fields. The deposition was done using e-beam 
evaporation of the corresponding films. The Ti02 tun¬ 
nel barriers were evaporated in 2 steps of 0.4 nm of Ti, 
followed by 15 minutes of oxidation in pure oxygen gas. 
The 65 nm thick Co contacts and the top gate were evap¬ 
orated in one step followed by an A1 capping layer. 

The sample configuration is shown in the inset of Fig. 

The bilayer graphene, in black, is encapsulated be¬ 
tween the top and bottom BN. The flake is fully encapsu¬ 
lated in the central region while both left and right sides 
are not encapsulated but only supported on a bottom 
BN. This configuration allows us to have ferromagnetic 
contacts at both sides of the sample while keeping the 
central region protected. This nonhomogeneous design 
also has consequences for the analysis of the data as we 
explain below. 

The top-gate together with the Si back-gate (Fig. 0(b) 
inset), allow us to simultaneously control the electric field 
E = ebgiVbg - V^°^)/2dbg - etgiVtg - v}g^)/2dtg and the 
carrier density n = eoCbgiVbg - vl,g^)/edbg + eoetg(Vtg - 

vl;g^)/edtg applied to the dual-gated region. Here e is the 
electron charge, eg is the vacuum permittivity, Cbg{tg) ~ 
3.9 the dielectric permittivity, dbg{tg) the thickness of the 

dielectric, Vbg(tg) the applied gate voltage and V^gltg) 
voltage at charge neutrality point of the back-gate (top- 
gate) respectively]!!]. 

We have characterized 3 devices (A, B and C) show¬ 
ing similar results at room temperature and 4 K. The 
results are shown in Table U From there we can see that 
Ts seems to depend on the length of the encapsulated 
region and, even though device G shows a higher spin 
diffusion coefficient (Dg) than device B, indicating bet¬ 
ter electronic quality, Tg of device B is more than 2 times 
longer than the one of device C. This nonstraightforward 
connection between electronic quality and Tg seems to 
be in agreement with the results shown in m for sin¬ 
gle layer graphene, while the length dependence can be 
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TABLE I. Spin parameters obtained at the gate combination 
giving the longest As for 3 devices with different length of the 
encapsulated regions. is the length of the encapsulated 
region and ci 2-3 is the separation between the inner contacts 


Dev. Lenc (pm) 

d2-3 (pm) 

T (K) 

rs(ns) Ds(m^/s) 

As (pm) 

A 13.2 

14.6 

300 

2.5 

0.07 

13 



4 

2.9 

0.2 

24 

B 8.5 

10.3 

300 

1.1 

0.03 

5.7 



4 

1.9 

0.05 

9.7 

C 6.2 

7.8 

300 

0.32 

0.04 

3.6 



4 

0.45 

0.07 

5.6 


explained by the effect of the invasive contacts and the 
lower quality of the nonencapsulated regions being re¬ 
duced increasing the contact separation [T3 [18]. As a 
consequence, we believe that even longer spin relaxation 
times can be achieved by making longer devices. 

From now on we will discuss the results obtained at 
4 K for device A. The contact resistances range between 
280 n and 2.7 kfl. These low values are a consequence 
of imperfect tunnel barriers and affect the spin transport 
measurements im. 

To characterize the charge transport properties of 
our device we carried out standard 4 probe measure¬ 
ments. We determine the resistances of the encapsulated 
and non-encapsulated regions by measuring the voltages 
(V 2 - 3 ) between contacts 2 and 3 and (^ 3 - 4 ) between 3 
and 4 while driving a current (/i-s) between 1 and 5 
(Fig. 0 (b)). In Fig. (a) we show the square re¬ 
sistance {Rsq) of the encapsulated region as a function 
of the back-gate voltage {Vbg) and the top-gate voltage 
(Vtg). The charge neutrality point appears as a line with 
a slope —CtgjCbg = —0.078 showing a resistance mini¬ 
mum at Vbg = —8 V, Vtg = —0.7 V. Taking into account 
that this point has zero carrier density and zero electric 
field we can estimate the doping at the top and bottom 
sides of the flake: « nfg wS.SxlO^® m~^. 

The resistance increases at both sides of the charge 
neutrality line, reaching up to 38 kfl at an electric field of 
0.69 V/nm. This is caused by the opening of a gap driven 
by the electric field m- One can also distinguish two 
Vtg independent features coming from the non-top-gated 
region between the inner contacts. One comes from the 
sides of the encapsulated regions that are non-top-gated 
and show a charge neutrality point around Vbg = —19 V. 
The other comes from the nonencapsulated regions and 
its square resistance is shown in Fig. acb). 

The mobility (^) obtained for this sample at Vbg = 
—8 V (zero electric field) is 16 m^/Vs at the electron side. 
The value is obtained using the formula Rgq = 1/ne^+ps 
where « 57 O is the n-independent resistivity compo¬ 
nent coming from the effect of the resistance of the non¬ 
top-gated regions between the inner contacts and short 
range scattering m- We have also confirmed that the 



FIG. 1. (a) Square resistance obtained between contects 2 
and 3 (in color scale) with respect to Vtg (y axis) and Vbg 
(x axis) (b) Square resistance of the non-encapsulated region 
measured between contacts 3 and 4. Inset: schematics of the 
device. 


resistance of the outer regions of the sample does not 
depend on Vtg. 

In Fig. Bb) we show the Vbg dependence of one of 
the outer regions’ resistance. As it can be seen from the 
graph, the charge neutrality point is below Vbg = — 50 V 
and falls outside our gate range. We attribute this to 
the contamination given by the polymers used during 
fabrication. 

To measure the spin transport properties of the en¬ 
capsulated region, we used the standard non-local ge¬ 
ometry [5D| sending a current (7i-2) between contacts 
1 and 2 and detecting the voltage (Va-a) between 3 
and 5. When applying an out-of-plane magnetic field 
the spins undergo Larmor precession. Measuring = 
V'i-bjl\-i while sweeping the magnetic field we ob¬ 
tained the so called Hanle precession curves. To elim¬ 
inate spin-independependent effects we have taken Hanle 
curves for parallel and antiparallel magnetic configura- 
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FIG. 2. Hanle precession curves obtained at Vbg = —50 V, 
Vtg = —5 V (top panel) and Vbg =50 V, Vtg =4.35 V (bottom 
panel) with the corresponding fitting curves and extracted 
spin parameters. 

tions of the inner contacts and substracted them i?„; = 
~ where jg nonlocal resis¬ 

tance in the parallel (antiparallel) magnetic configuration 
[HI- The magnetizations of the contacts are tuned ap¬ 
plying an in-plane magnetic field. 

In Fig. we show two Hanle curves taken at Vbg = 
50 V, Vtg = 4.35 V and Vtg = -50 V, Vtg = -5 V, 
corresponding to the top right and bottom left corners 
in the color plot of Fig[^ (a). The spin relaxation time 
and diffusion coefficients are extracted from these curves 
by fitting them with the solution of the Bloch equations 
[Il|20] including a small offset [2T| . 

The spin signal at Vtg = 50 V, Vtg = 4.35 V (Fig. [^top 
panel) is 10 times higher than the one at Vbg = —50 V, 
Vtg = 5V (Fig. bottom panel). This is most likely due 
to the low resistance of our contacts. When the contact 
resistance is in the order of the spin resistance of the 
channel {R\ = RsqXs/W, here W is the width of the 
sample), part of the injected spin accumulation relaxes 
back to the contacts instead of diffusing in the channel, 
reducing the effective injection efficiency [HI HI] from 
12% to 2%. This effect is ruled by the resistance of the 
outer regions (where the contacts are placed) and it is 
Vtg independent. 

In Fig. (a) we show the spin and charge diffu¬ 
sion coefficients for 3 different backgate voltages in dots 
and solid lines respectively. The charge diffusion coef¬ 
ficients (Dc) were obtained using the Einstein relation 
l/Dg = e^Rsqv{Ep), where v{Ef) is the density of states 
at the Fermi energy and e the electron charge. Rsq was 
taken from Fig. (a) and corrected by substracting the 
resistance of the nonencapsulated regions between the in¬ 
ner contacts. 

At Vbg = — 50 V one can observe a substantial dif¬ 
ference between Dg and Dg- We attribute this to the 
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FIG. 3. Spin transport parameters obtained by fitting the 
Hanle curves at 4 K as a function of Vtg for Vbg = -50,-10 
and 50 V. (a) Spin diffnsion coefficient (dots) compared with 
the charge diffusion coefficient (solid lines), (b) spin relaxation 
time and, (c) the relaxation length. The lines connecting the 
spin parameters are a guide to the eye. 


fact that the gate induced doping of the encapsulated 
and nonencapsulated regions have different signs, creat¬ 
ing pn junctions of unknown widths at the boundaries. 
Since these boundaries are between the inner contacts, 
the measured square resistances are no longer character¬ 
istic of the channel itself but of the junctions. This affects 
the determination of Dc- 

At Vbg = 50 V both encapsulated and non¬ 

encapsulated regions are electron-doped and Dc shows 
better agreement with respect to Hg, supporting the va¬ 
lidity of the parameters obtained from the Hanle mea¬ 
surements. The minor discrepancy in this case can be at¬ 
tributed to the small resistance of the encapsulated non¬ 
top-gated region that was not substracted from the calcu¬ 
lation of Rsq. Dc and Dg reach values above 0.5 m^/s, 5 
times higher than the best achieved for BN encapsulated 
single layer graphene spin valve devices |12j . 
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At Vbg = —10 V (approximately zero electric field) we 
see that close to the charge neutrality point (Vtg « 0 V) 
there is a better agreement between Dc and Dg than 
at high carrier densities. This can be explained tak¬ 
ing into account that close to the charge neutrality point 
the square resistance of the encapsulated region is high 
enough to dominate the measurement of Rsq, but at high 
carrier densities, the square resistance of the encapsu¬ 
lated region is small and the contributions of the non¬ 
top-gated regions become relevant. 

In Fig. 1^ (b) there is a strong dependence of the spin 
relaxation time on Vbg. At Vbg = — 50 V, the relaxation 
time reaches 3 ns while at 50 V a maximum of Tg =310 ps 
is obtained. This reduction of a factor 10 in Tg is in 
agreement with the results in |12j and can be explained 
as an effect of the change in R\ of the nonencapsulated 
regions. As the spin resistance of these regions increases, 
their influence on the spin relaxation is reduced. This 
effect occurs because the spins relax predominantly at 
the regions with the lower Rx. 

The opposite effect occurs when opening a gap in the 
encapsulated region and its square resistance increases 
with respect to the one of the non-encapsulated part. 
Since Tg is longer in the encapsulated region than in the 
outer ones, Rx gets orders of magnitude higher than the 
one of the outer part. As a consequence, the spin re¬ 
laxation is dominated by the non-encapsulated regions 
and the amplitude of the spin signal vanishes. This ef¬ 
fect explains why we could not measure spin signals at 
Vbg = — 50 V and positive Vtg. 

In Fig. (c) we show the spin relaxation lengths cal¬ 
culated using the formula As = x/DgTg. Apart from its 
carrier density dependence, we can see that As goes up 
to 24 pm, the highest value achieved up to date by fit¬ 
ting Hanle measurements in nonlocal geometry. Note 
that, even though spin relaxation lengths up to 280 pm 
were estimated from local 2 probe measurements at 4 K 
for epitaxial graphene on a SiC substrate [53], no Hanle 
measurements were done to support these values. 

Since the spins probe the whole device (inner and outer 
regions), we have to account for the nonhomogeneity of 
our sample to explain our results |12j . For this reason, 
we use the same model as in [12] and [24], where we 
solve the Bloch equations for a nonhomogeneous system 
consisting of a central region sandwiched by 2 regions 
as shown in the inset of Fig. [^ We set the spin and 
charge transport parameters (r^, Dg and Rgq) for the 
3 regions assuming that the outer regions are identical. 
After simulating the corresponding Hanle curves, we fit 
them using a homogeneous model as we have done with 
our experimental data. From this procedure, we obtain 
the effective relaxation time of the system (re//). 

The effect of invasive contacts is taken into account us¬ 
ing the spin transport parameters obtained from Hanle 
measurements carried out at the non-encapsulated re¬ 
gions. Since the contact separation in these regions is be- 
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FIG. 4. Effective spin relaxation time in the system as a 
function of Tenc for 3 different values of r„on at Vbg = —50 V 
and Vtg = —5 V (a) and Vbg = 50 V and Vtg = 4.35 V (b), 
the dashed line is the experimental value of Te//, taken from 
Fig. ^ (top panel). The inset shows a cartoon device of the 
simulated system. 


tween 1 and 2 pm, the extracted parameters are strongly 
affected by the low contact resistances m 

In Fig. I^we plot Teff as a function of the spin relax¬ 
ation time in the encapsulated region (Tenc) for 3 dif¬ 
ferent values for the spin relaxation time in the non- 
encapsulated region (r„o„). The resistance values used 
for the central region are the ones used to calculate Dc 
in Fig. (a). For the diffusion coefficient in the en¬ 
capsulated region {Dgnc) we have used the values of Dg 
extracted from the experiments at the encapsulated re¬ 
gion. This is justified since Dg is not affected by the 
outer regions [12] [24]. The square resistance of the non- 
encapsulated region is taken from Fig. (b) and the dif¬ 
fusion coefficient {Dnon) is taken from the experimental 
Hanle curves obtained at the outer region. 

In Fig. (a), where Vbg = —50 V and Vtg = —5 V, 
the maximum Tg// obtained from the simulations reaches 
1.8 ns for Tenc = 100 ns and r„o„ = 300 ps. This value 
is still below the 2.9 ns obtained from the fittings of the 
encapsulated data. This discrepancy can be explained 
taking into account that, as discussed above, there are 
pn junctions in the boundaries between encapsulated and 
non-encapsulated regions and this is not taken into ac- 
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count in our simulations. 

The simulations at Vbg = 50 V and Vtg = 4.35 V are 
shown in Fig. (b). The dashed line at 310 ps cor¬ 
responds to the value obtained from the fittings of the 
experimental results at the encapsulated region. The in¬ 
tersections between the simulated curves and the dashed 
line give us the possible values of Tenc- From the fit¬ 
tings of the Hanle curves measured at the outside re¬ 
gions, we obtained Tnon ~100 ps. As a consequence, 
from our simulations, Tenc ~ 1 ns. Using this relaxation 
time and Dene = 0.52 m^/s, we obtain a spin relaxation 
length of 22 pm, close to the 24 pm at Vbg = — 50 V and 

Vtg = -5Y.^ 

In conclusion, we have characterized 3 boron nitride 
encapsulated bilayer graphene devices with 13.2, 8.5 and 
6.2 pm long encapsulated regions showing consistent be¬ 
havior. 

Our results show that the measured Tg depends on the 
length of the encapsulated region in agreement with |12j . 

The results obtained for the longer device show 
unprecedented large spin diffusion coefficients up to 
0.52 m^/s at 4 K, 5 times higher than the best achieved 
for single layer graphene using the same geometry m- 
As a consequence, the spin relaxation length rises up to 
13 pm at room temperature and 24 pm at 4 K. 

Our simulations using a three regions model show that 
the measured spin relaxation times of 2.5 ns at room tem¬ 
perature and 2.9 ns at 4 K are most likely limited by the 
outer regions, suggesting that it is possible to transport 
spin information over even larger distances in the used 
geometry by increasing the length of the encapsulated 
region. According to this result, higher spin relaxation 
times can be achieved by making longer devices [15] . 
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